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In this work, the many-spin interactions taking place in 
Mni2 large-spin clusters are extensively studied using the 8- 
spin model Hamiltonian, for which we determine the pos- 
sible parameters based on experimental data. Account of 
the many-spin excitations satisfactorily explains positions of 
the neutron scattering peaks, results of EPR measurements 
and the temperature dependence of magnetic susceptibility. 
In particular, strong Dzyaloshinsky-Morya interactions are 
found to be important for description of neutron scattering 
data. The role of these interactions for the relaxation of the 
magnetization is qualitatively discussed. 
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INTRODUCTION 

In the last years, a new kind of magnetic compounds, 
the magnetic molecules, has bees drawing the attention 
of physicists as well as chemistsEl. Such molecules each 
contain a large number (typically, 10 to 20) of para- 
magnetic ions (such as Mn, Fe or Cu) coupled by ex- 
change interactions. Each molecule, therefore, presents 
a mesoscopic system that is neither totally microscopic, 
nor totally macroscopic, but where micro- and macro- 
scopic behavior coexist. These materials are promising 
for various practical applications!! On the other hand, 
the coexistence of quantum and classical behavior in the 
clusters makes them very suitable objects fcajj-study of 
macroscopic quantum effects in spin systemsBu. These 
studies, clarifying many problems of quantum theory of 
measurements^, are also important for development of a 
physical basis for practical implementation of powerful 
algorithms of quantum computations, quantum cryptog- 
raphy and quantum searchingQ. 

Particularly, the Mni 2 Oi2(CH 3 COO)i 6 (H20)4 mole- 
cules (below referred to as Mni2) recently became a sub- 
ject of great interest. Each moleculeQiS contains a cluster 
of twelve manganese ions surrounded by acetate radicals 
and water molecules. The ground state of the clusters 
corresponds to a large total spin S — 10. The clusters 
possess a strong easy-axis anisotropy: the zero-field split- 
ting between the states with S z = ±10 and S z = ±9 
(where S z is the value of z-projection of the total cluster 
spin) is 14.4 K. Being stacked into a crystal, the molecules 
form a tetragonal lattice; in so doing the magnetic inter- 
actions between different clusters are very small (of order 



of 10~ 2 T). Thus, the crystal consisting of these molecules 
can be considered as an assembly of ideal noninteracting 
superparamagnetic entities, each being identical to the 
others. 

These clusters have been successfully used for the 
study of mesoscopic quantum effects. In particular, res- 
onant magnetization tunneling has,r-byeen unambigously 
registered in experiments xtp Mni2ull3. Moreover, there 
are experimental resultsEj'Ej supporting the hypothesis of 
"ground state-to-ground state" tunneling in Mni 2 below 
2 K (for more detailed discussion, see Sec. VI). 

However, the progress in understanding the physical 
properties of Mni2 is greatly hampered by the lack of an 
adequate description of these clusters. Indeed, the de- 
scription of Mni2 as a single spin S = 10 entity has been 
the starting point in most works devoted to this sub- 
ject. We know of only a few theoretical attempts r ±p-.ac | - 
count for the internal spin structure of the clusterO'oEj, 
but even in these the relativistic anisotropic interactions 
have not been— taken into account. In view of recent 
experimental30Ea showing that the single-spin model 
is seriously deficient, it is worthwhile reconsidering the 
many-spin aspects of Mni2. 

In this paper we focus on the many-spin interactions 
in Mni2 clusters. We account for not only isotropic 
exchange interactions between ions in the cluster, but 
also various anisotropic interactions possibly present in 
Mni2. Based on the results, we propose a spin Hamil- 
tonian for these clusters. We show that this Hamilto- 
nian can reproduce satisfactorily most recent experimen- 
tal results, such as positions of neutron scattering peaks, 
high-frequency EPR data and the experimental depen- 
dence of the magnetic susceptibility on temperature. We 
note that the account of anisotropic interactions, espe- 
cially the Dzyaloshinsky-Morya interaction (which has 
been missing up to now) , is crucial for a detailed descrip- 
tion of the experimental data. 

The paper is organized as follows. In Sec. | we de- 
scribe the basic model of Mni2 used in this work and 
establish roughly its domain of validity. In Sec. || we 
deri ve a nd discuss the spin Hamiltonian for this model. 
Sec. HI is devoted to discussion of relevant experimental 
data. In Sec. IV the numerical procedure used for calcula- 
tions is discussed and the possible parameters of the spin 
Hamiltonian are presented. Comparison with experimen- 
tal data is made. The results obtained are analyzed quali- 
tatively and discussed in Sec. |y|, where the interpretation 
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of the neutron scattering data is presented. In Sec. VI we 
qualitatively discuss the relation between Dzyaloshinsky- 
Morya interactions and mag netic relaxation in M2112. A 
summary is provided in Sec. VII . 



I. THE DIMERIZED 8-SPIN MODEL OF Mn u 

The cluster Mni2, schematically shown in Fig. |l], con- 
sists of eight Mn 3+ ions having the spin 2 and four Mn 4+ 
ions having the spin 3/2. The ions are coupled by ex- 
change interactions, indicated in Fig. [| by different lines 
connecting the ions. The values of the exchange inte- 
grals are not known, but estimates are given in Ref. @: 
Ji = -150 cm" 1 (AFM exchange), J 2 = J3 = -60 cm^ 
and J 4 < 30 cm -1 . These values are rough, but describe 
correctly the scale-of-exchange interactions in Mni2- Re- 
cent experimentsOEZHia show that the excitations with 
spin values S < 10 are rather close to the ground state: 
the distance is 40-60 K (values differ in different reports) . 
This is less than the energy of some states with the spin 
S = 10 (namely, the states S z — 0, ±1, ±2, ±3), i.e. the 
lower states of the manifold 5 = 9 are lower than the 
higher states of the manifold S = 10. Thus, an adequate 
description of Mni2 should account for the excitations 
with S < 10; i.e. the cluster should be considered as a 
many-spin system. 

The total number of spin states in Mni2 is large even 
for modern computers. But we can employ the fact that 
the exchange antiferromagnetic interactions J\ (see Fig. 
|l|) are much larger than all the othersa, so corresponding 
pairs of ions Mn 3+ and Mn 4+ form dimers with the to- 
tal spin s = 1/2 (one of this pairs is designated in Fig. 
[l], it includes ions C and D). This model has already 
been successfully used for description of spin states of 
the clusteiSpl Its validity is proven by megagauss-field 
experimental^: the states of dimers with the spin s higher 
than 1/2 (excitations of dimers) come into play when the 
external magnetic field is about 400 T, i.e. the excitations 
of dimers have energy about 370 cm . Analogously, the 
dependence of the inagnetic susceptibility of the cluster 
versus temperatureO'EH shows that the dimer excitations 
contribute when temperature becomes as high as 150-200 
K. 

Based on these data, we can analyze the domain of 
validity of the "dimerized" model. To do this, we note 
that the exchange interactions J 2 , J3 and J4 mix the 
ground state of a dimer with the dimer excitations, and 
the approximation of spin- 1/2 dimers corresponds to the 
zcroth order perturbation theory with 1/ J\ as an expan- 
sion parameter (similar approach has been used in Ref. 
0). 

To clarify this point, let us consider the level a having, 
to the zeroth order, the energy E a with respect to the 
ground state. Let us denote the distance between the 
ground state and the excitation of dimer as E ex ~ 370 
cm . The first-order correction to the energy of the 



level a is of the order of J' 2 / (E ex — E a ), where J' is the 
magnitude of exchange interactions between dimers and 
nondimerized spins (see below). Thus, accounting for the 
first-order corrections, the distance between the ground 
state and the level a becomes 

E' a = E a + C a J' 2 [l/{E ex - E a ) - 1/E„)] , 

where C a is a factor of order of unity, depending on the 
specific level a. As will be shown below, J' is of order of 
70 cm -1 ; so the first-order correction for the levels with 
energies about 70 cm -1 is already considerable, of order 
of 4 cm -1 . This estimate, though being rough, gives the 
correct order of magnitude of the error introduced by the 
dimerized 8-spin model. 

Moreover, this error restricts the region of tempera- 
tures where the dimerized model can be successfully ap- 
plied. E.g., as our calculations show, to obtain the correct 
value of the magnetic susceptibility \ at the tempera- 
ture T, we need to account for the levels with energies 
about 4-5 kT. Obviously, the error in positions of these 
levels will introduce corresponding error in the depen- 
dence x(T). Its analytical evaluation is difficult, and the 
comparison of the results of calculations with the exper- 
imental data, performed in Sec. IV is the better way to 
understand the temperature domain of validity of the 
dimerized model. As our results show, the dimerized 
model gives reasonable results for temperatures lower 
than about 50 K. 

Recalling that the temperatures below 30 K are of most 
interest, we conclude that the dimerized model is satis- 
factory for present needs of experimentalists. 



II. THE SPIN HAMILTONIAN OF Mn 12 

Thus, we consider the Mni 2 cluster as consisting of four 
" small" dimer spins s = 1/2 and four " large" spins 5 = 2 
(corresponding to the four non-dimerized ions Mn 3+ ), 
coupled by exchange interactions (see Fig. ||). Moreover, 
we have to account for the anisotropic relativistic inter- 
actions in the cluster, so the Hamiltonian of the system 
can be written as: 



n 



-j 



J' ^ SfcS; + H Te \, 

(k,l) 



(1) 



where Sj are the spin operators of small dimer spins 
s = 1/2, S/ are spin operators of large spins S = 2, 
and H rc \ denotes the part of the Hamiltonian describing 
relativistic interactions in the cluster. Summation in (|l|) 
is over pairs of spins coupled by exchange interactions. 
In the first term of the Hamiltonian we took into account 
that each small dimer spin is coupled with all the other 
small spins, so 2J2 s i s j = E s i) U P t° an insufficient 
additive constant. 

To zeroth order in J%, the exchange integrals of the 
dimerized models are connected with the initial exchange 
parameters J2, J3 and J4 as follows: 
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J 



-J2/2, 



J' 



- J 3 + 2J 4 



(2) 



Since the values of J2, J3 and J4 are not known, the 
parameters J and J' are to be determined from experi- 
mental data (see Sec. IV). 

Furthermore, different types of relativistic anisotropic 
magnetic interactions possibly present in Mni2 clusters 
should be included in the Hamiltonian. A large easy-axis 
anisotropy in the cluster is one of most important fea- 
tures to be taken into account. Generally, this anisotropy 
arises due to the single-site anisotropy of large spins 
(spins of Mn 3+ ions) and various kinds of anisotropic ex- 
change. We performed calculations for three basic types 
of easy-axis anisotropy in the cluster: 
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(3a) 
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(3b) 
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-Jzz ^ sfSj, 

(id) 


(3c) 



where summations in (3b) and ( |3c| ) are over exchange- 
coupled pairs of spins. Anisotropy parameters (K z , J zz 
or Jzz) have been chosen to give a correct value of the 
zero-field splitting between the states S z = ±10 and S z = 
±9 (14.4 K). All three types of anisotropy give rather 
close energies of low- lying excitations (of energy less than 
40 K), but higher excitations are reproduced best if the 
anisotropy is assumed to be of single-site type ([3a|) , so we 
can conclude that the easy-axis anisotropy is primarily 
of single-site type. This result agrees with the conclusion 
drawn in Ref. |2l]. We will consider only this kind of 
anisotropy. 

Another potentially important sort of relativistic inter- 
action is an in-plane anisotropy of large spins, i.e. H Te \ 
can include a contribution of the form: 



(S v 2 ) 



(4) 



where the presence of fourth-order symmetry axis in the 
cluster is directly taken into account. The small spins 
s = 1/2 are excluded since (a x ) 2 = (<x y ) 2 = {<r z ) 2 = 1 
for Pauli matrices a x , a y and er z ; and only spins of 
non-dimerized Mn 3+ ions give a nontrivial contribution. 
These ions are surrounded by eight oxygen ions form- 
ing a distorted octahedron. The axes of oxygen octahe- 
dra are significantly tilted from the c-axis of the cluster, 
therefore, this term can be relatively large, even compa- 
rable to the easy- axis anisotropy. But, surprisingly, our 
results show that this kind of interaction gives negligi- 
ble effect, except for trivial renormalization of the easy- 
axis anisotropy constant K z in (3a). If we account for 
this renormalization, the positions and the wave func- 
tions of excited levels remain almost unaffected even for 
K\ = 3K Z (i.e., for the in-plane anisotropy three times 



larger than the easy-axis one). Thus, this kind of inter- 
action can be excluded from further considerations. 

Another important interaction is Dzyaloshinsky-Morya 
(DM) antisymmetric exchange. To our knowledge, the 
possible presence of DM-intcractions in Mni2 was first 
suggested in Ref. [l8[ but little attention has been paid 
until now. Our results show that these interactions are, 
indeed, very important and have rather large magnitude. 

A pair of ions coupled by DM-interaction is described 
by the Hamiltonian 



DM 



D - [Si x S 2 ], 



(5) 



ancLthc magnitude of the DM- vector D can be estimated 
asE3 D ~ \A, where A is the isotropic (nonrelativis- 
tic) exchange coupling between ions and A is the spin- 
orbit coupling constant (which is rather small for transi- 
tion ions). For comparison, the magnitude of easy-axis 
anisotropy is estimated as K z ~ X 2 A, i.e. is of next or- 
der of smallness in comparison with D. Thus, the DM- 
interactions in Mni2 can be expected to be important. 

In the 8-spin model of the cluster there are DM- 
interactions of two kinds: 



(id) 

ff DM = EE D i i '[ s < xs i] 



(6a) 
(6b) 



Summation in (|6a| ) is over exchange-coupled pairs of 
spins; summation in ( |6~b| ) is over all pairs of dimer spins, 
since all dimer spins interact with each other. We studied 
both kinds of DM-interaction and found that the second 
kind, i.e. H^ M involving small spins can be neglected. 
Therefore, we can neglect the interactions of the type 

The crystal field in Mni2, governing the DM- 
interactions, possesses certain symmetry elements, thus 
imposing restrictions on the values of D JJ . It is reason- 
able (and rather standards) to assume that the crystal 
field is determined mainly by the oxygen octahedra sur- 
rounding manganese ions in the cluster; so the symmetry 
of the crystal field is governed by the mutual arrangement 
of the oxygen octahedra. The following two symmetry el- 
ements are of interest for us n The first one is the fourth- 
order rotary-reflection axisErB parallel to the c-axis of the 
cluster. This symmetry is obviously preserved in the 8- 
spin model of the cluster, so the two DM-vectors D 1,5 
and D 1 - 8 (see Fig. |) define all the other D*'*- 7 '. The other 
element of symmetry is the mirror plane p parallel to 
the z-axis passing through the ions CandD (see Fig.|). 
The oxygen octahedra surrounding the ions A and B (see. 
Fig. 0) are invariant with a good degree of accuracyEHl 
with respect to reflection in the plane p (inspection of 
the structure data supplied in Refs. ^ shows this); this 
symmetry is also preserved in the 8-spin model. Thus, 
the vector D 1,8 (Fig. 0) defines all the other DM-vectors 
in the Hamiltonian (pal): 
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Dl* = -D\* = D 2 / = -D 2 / = -D 3 f = D 3 / (7a) 

= -Dtf = Df, 
D 1 / = D 1 / = -D 2 / = -D 2 / = -D 3 / = -Df (7b) 

= D 4 / = D 4 /, 
D 1 / = -Dl> 5 = D 2 / = -D 2 / = Df 6 = -D 3 / (7c) 

= D iJ = -Di>*. 

Obviously, any other vector D 4J can be taken as a ba- 
sis instead of D 1 ' 8 . No other symmetry elements allow 
for further reduction, so DM-interactions in Mni2 are 
described using three parameters: D x ' 8 , D y ' & and D\' & . 
Below, these parameters are denoted simply as D x , D y 
and D z . 

As our results show, in the DM-Hamiltonian ( |6a|) the 
terms proportional to D y produces negligible matrix el- 
ements (a few percent in comparison with other terms). 
It occurs due to symmetry reasons: the inspection of 
the relations ([[]) shows that the components D x and D z 
transform antisymmetrically with respect to reflection in 
the plane p, but the component D y transforms symmet- 
rically. The matrix elements of the terms proportional 
to D y nearly cancel each other, leading to negligible ma- 
trix elements. Therefore, these terms are excluded from 
consideration and we set D y = with negligible error. 

Finally, having studied all the interactions described 
above, we can write down the Hamiltonian of the cluster 
in the following form: 

n = - s *) 2 - J ' E s ^ - K * E ( 5 ?) a ( 8 ) 

i (k,l) i=l 

+ £D«.[s i xS i ], 

(id) 

where the DM-vectors D lJ obey the relations (^) with 
the parameter D y ' s = D y = 0. 



III. REVIEW OF RELEVANT EXPERIMENTAL 
RESULTS 

At present, data of various experiments on magnetic 
molecules Mn^Ac are available, including the tempera- 
ture dependenge,p£ the effective magnetic moment oi-.the 
cluster /i c fj(T)BEni3, the results of EP-B. experimentsEJS, 
dynamic susceptibility measurements^, inelastic neutron 
scattering datat-l and specific heat datatll. Unfortu- 
nately, only few of these data can be used for determining 
the parameters of the 8-spin Hamiltonian for Mni2 clus- 
ters, f— | 

Recent high-frequency EPR experiments^!] refined the 
description of the easy-axis anisotropy of the cluster and 
showed that the anisotropy Hamiltonian in the single- 
spin model can be approximated as follows: 

H = aS 2 + (3St + i(Sl + Si), (9) 



a = -0.56K, (3 = —11.08 • 1CT 4 K, 
7 = 2.88-l(T 5 K, 

where S z , S+ and S- denote the operators of the total 
spin of the cluster. It means, in terms of a many-spin 
approach, that the energies of the low-lying levels with 
spin S = 10 obey Eq. (|9|). It is worth noting that the 
derived values of quartic corrections P and 7 are rather 
large and, as our calculations show (see below), seem to 
be poorly explained using the single-spin description of 
Mni2, i.e. when accounting only for the states belonging 
to the S = 10 manifold. Our results show that the excited 
levels with S < 10 are necessary to give reasonable values 
for the quartic corrections. 

Another set of results, very useful for elucidating the 
many-spin interactions in Mn^Ac, is the neutron scat- 
tering results supplied in Ref. [l^. The experiments have 
been performed at very low temperatures (mostly, 1.5 K 
to 2.5 K), where only the lowest levels S z = ±10 are pop- 
ulated. Since the selection rule for neutron scattering is 
SS Z = 0, ±1, only the levels with S z = ±9 can give rise 
to scattering peaks (the levels with S z = ±11 have too 
large energies^ and can be excluded). 

Results of these experiments can be summarized as fol- 
lows. A prominent peak of spin origin at about 0.3 THz 
has been detected and attributed to the transitions to the 
levels with S = 10, S z = ±9, in excellent agreement with 
all previous data (0.3 THz corresponds to about 14.4 K). 
At higher energies, two sets of peaks have been detected 
around 1.2 THz and 2.0 THz. The fitting proposed in 
Ref. [l6] gives two peaks in the first set (at energies about 
57 K and 66 K) and three peaks in the second set (at 
energies 90 K, 96 K and 105 K); but authors indicate 
clearly that possibly more peaks are present (most likely, 
three peaks in the first set and four or five in the second). 
Another interesting detail of the neutron scattering spec- 
tra is a very broad mode situated at about 0.2 THz; this 
mode disappears when the temperature is less than about 
2 K. 

The authors have not managed to interpret these fea- 
tures, except for the peak at 0.3 THz. They pointed out 
that there is particular difficulty in interpretation of the 
peaks at 1.2-1.3 THz: the model they used for suscepti- 
bility fitting gives two degenerate levels S = 9 at about 
33 K, an obvious contradiction with the neutron scat- 
tering spectrum. We show below that the 8-spin model 
developed here can overcome these difficulties and gives 
correct positions for neutrons peaks at 1.2 THz along 
with a correct description of the susceptibility data. 

Thus, we found the following experimental results to 
be relevant for the purpose of a quantitative description 
of the Mni2 clusters. The distance between the ground 
state and the first excited level(s) is 14.4 K. The energies 
(the anisotropy splittings) of the low-lying levels, belong- 
ing to the S = 10 manifold, obey formula (|9|). There are 
two or three neutron peaks around 60-70 K, two of them 
are situated at 57 K and 66 K. Also, there are up to five 
peaks around 100 K, three of them are at 90 K, 96 K and 
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105 K. The temperature dependence of the susceptibility 
(or, equivalently, the dependence fj, c s(T)) has the form 
displayed in Refs. g, |, |l| and Fig. |. 

The other experimental results, though providing im- 
portant information about Mni2, are much less suitable 
for our purposes (to a large extent, because different, a 
priori equally probable, interpretations are possible). 



IV. NUMERICAL CALCULATIONS AND 
PARAMETERS OF THE 8-SPIN MODEL. 
COMPARISON WITH EXPERIMENT 

Having derived the spin Hamiltonian for the 8-spin 
model of Mn 12 , we attempted to extract its parameters 
from the relevant experimental data. 

We used the following two-step numerical scheme. At 
the first step, the relativistic term H le \ has been ne- 
glected resulting in an isotropic exchange Hamiltonian. 
The eigenstates of this Hamiltonian are degenerate with 
respect to S z . Thus, it is sufficient to take into account 
only the states with S z — 0, so the exchange Hamiltonian 
(represented by a matrix 1286x1286) has been diagonal- 
ized within the subspace spanned by these states. Then, 
at the second step, the relativistic anisotropic interac- 
tions have been taken into account. Among the states 
obtained at the first step (having S z = 0) , we retain only 
those with the energy less than E cut (a sufficiently large 
value for this parameter has been chosen) and generate 
the corresponding states with different S z (basis states). 
Then the complete Hamiltonian (||) has been diagonal- 
ized within the subspace spanned by the generated basis 
states. Calculations with different values of _E cut have 
been performed to assure that the positions of lower lev- 
els are obtained with desired accuracy. Typical values of 
-Scut were about 250 K: the levels with higher energies 
are not worth including due to the limited accuracy of 
the 8-model itself (see Sec. |). 

Based on the procedure described above, the fitting 
of relevant experimental data (Sec. HI) has been made 
and the possible parameters of the 8-spin Hamiltonian 
determined. 

Neutron scattering data are of primary interest for us. 
We focus our attention on the positions of the neutron 
peaks, since the amplitudes depend strongly on details 
of the experiments. We first assume an ideal experi- 
ment, where the resolution of the setup is infinite and 
the neutrons with all possible scattering vectors are de- 
tected (i.e., the detector has infinite aperture). In this 
case, at zero temperature, tha-cross-section of neutron 
scattering at the energy E isc3c3 



\SW)W\S b n \0)6m)-E), 



(7(E) = [ d 3 q AF 2 (q) V (S a , b - q a q b /q 2 ) (10) 



where A is a constant, F(q) is the form-factor of man- 
ganese ions, q is a scattering vector, n, m enumerate 
different ions, and a, b refer to the Cartesian coordinates 
(x, y and z). The integration is performed over all vectors 
q. E(tp) denotes the energy of the state tp, |0) denotes 
the ground state, which is the only one populated at zero 
temperature. The transitions with AS Z = 0, +1 can be 
neglected since there is only one state S z = 10 and no 
states S z = 11. In this case, the total cross-section ( |l0| ) 
is proportional to the quantity 



(ii) 



where the state ifi has the energy E(ip) = E with re- 
spect to the groundstate; i.e., the state ip is the final 
state of the neutron scattering process and gives rise to 
a neutron peak at the energy E(ip) of the amplitude pro- 
portional to V. The summation in (11) is performed 



J^exp [«q(r m - r„) 



over all basis levels having S z = 9; these levels are de- 
noted as </>' 9 '. Equation ( pi] ) expresses the simple fact 
that only the transitions with AS Z = — 1 are allowed in 
the neutron scattering process, since the transitions with 
AS Z — 0, +1 are absent. Below, the quantity V is re- 
ferred to as a normalized cross-section for the level if). 
Our results show that V discriminates easily the eigen- 
states which can give rise to noticeable neutron scattering 
peaks. 

Furthermore, the values of parameters a and (3 de- 
scribing the easy-axis anisotropy in Eq. (^J), have been 
taken into account in determination of the cluster pa- 
rameters. The energies of the five lowest levels, having 
spin S = 10, have been approximated by a fourth-order 
polynomial, following Eq. (^), and the coefficients a and 
(3 have been extracted and compared to the experimental 
data. 

As a result of calculations, the following three sets of 
the cluster parameters have been found to provide the 
best fitting of experimental data: 

Set A: J = 0, J' = 105 K, K z = 5.69 K, D z = -1.2 K, 
D x = 25 K; 

Set B: J = 23.8 K, J' = 79.2 K, K z = 5.72 K, D z = 10 
K, D x = 22 K; 

Set C: J = 41.4 K, J' = 69 K, K z = 5.75 K, D z = 10 K, 
D x = 20 K. 

The positions of neutron peaks calculated for these sets 
of parameters are presented in Fig. |[ The graphs show 
the dependence of normalized cross-section vs. the level 
energy. It is seen from these figures, that the normal- 
ized cross-section is extremely small (less than 10 -2 ) for 
most of levels, and only few states can give rise to no- 
ticeable neutron scattering peaks. Moreover, to facilitate 
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the analysis of the data for reader, the positions of neu- 
tron peaks are listed in the Table 0. The values of the 
easy-axis anisotropy parameters a and (3 are listed in the 
Table || 

As the results show, each of the parameter sets repro- 
duces reasonably well its own portion of the experimen- 
tal results. All the sets give reasonably good positions of 
the low-energy neutron peaks at 0.3 THz (14.4 K), 1.19 
THz (57 K) and 1.38 THz (66 K). The parameter set A 
also gives the values of anisotropy parameters a and /3, 
rather close to the experimental ones, but the neutron 
peaks corresponding to higher energies (around 2 THz) 
are reproduced poorly. The parameter sets B and C give 
correctly only the order of magnitude of a and (3, but re- 
produce better the positions of the high-energy neutron 
peaks. 

Finally, the temperature dependence of the effective 
magnetic moment /z e ff of the cluster has been calculated 
for all three sets of parameters using the formula 



Hes(T) = y/3 X (T)-kT, 



(12) 



where x is the susceptibility of the cluster, k is Boltz- 
mann's constant and T is the temperature. The suscep- 
tibility x(T) has been calculated in a way reproducing 
the experimental procedure. The Zccman term, describ- 
ing the effect of an external field has been introduced into 
the Hamiltonian (|§|) . The field magnitude H = 1 mT has 
been chosen following Ref . [l6| The resulting Hamiltonian 
has been diagonalized, and the component of the cluster 
spin 6>h along the field has been calculated by means of 
quantum-statistical averaging over the Gibbs canonical 
ensemble. This routine has been repeated several times 
for different orientations of the field, and the obtained 
values of iSh have been averaged. It corresponds to a 
powder sample measurements, when the crystallytes are 
randomly oriented with respect to the field. The suscep- 
tibility x, following a standard experimental procedure, 
has been calculated as a ratio of the resulting average 
cluster spin to the field magnitude H. Finally, the value 
yix e ff has been obtained applying Eq. (|l2|). 

The curves [i c s calculated for the three sets of clus- 
ter parameters, are presented in Fig. |^ along with the 
experimental data. All the sets give almost coinciding 
curves, and below 50 K the agreement with experiment 
is good. The region of temperatures higher than 50 K can 
not be reproduced satisfactorily: as our test calculations 
showed, to obtain the correct value of effective moment 
/x e ff at the temperature T, we need to account for the 
levels with energies about 4-5 kT. When calculating the 
curves presented, only the levels with energies less than 
250 K have been taken into account, thus restricting the 
correctly described temperature region. 



V. QUALITATIVE ANALYSIS OF THE RESULTS 
AND INTERPRETATION OF EXPERIMENTAL 
DATA 

At present, having rather limited number of the rele- 
vant experimental data, it is hard to distinguish between 
the parameter sets A, B and C. The easy-axis anisotropy 
parameters a and (3 are obtained with good precision in 
EPR-yxpcriments, but the magnetization measurement 
datalia suggest other values for these parameters; so, com- 
parison of the experimental values of a and (3 with our 
results can not serve as a definitive basis for judgement. 
Also, the quality of the description of the high-energy 
neutron peaks can not be decisive, since the disagreement 
can be attributed to the limited accuracy of the dimpt- 
ized 8-spin model itself. Megagauss-field experimentsEj, 
along with careful measurements of the low-energy peaks 
(around 1.2 THz) and fitting of their amplitudes seems 
to be a promising strategy for future investigations. 

Nevertheless, the results already obtained provide new 
important information about the role of many-spin in- 
teractions in Mni2 clusters. In this section we focus our 
attention on the qualitative consideration of the results 
obtained and discuss the interpretation of the experimen- 
tal data. 

First, it is worthwhile to note that the consideration 
of states with spin <S less than 10 leads to rather large 
quartic corrections to the energy of easy-axis anisotropy. 
If these excited states are not taken into account, i.e. if 
only the states with S = 10 are included, the value of (3 
is of order of 10~ 5 K. 

Another important fact is the large magnitude of the 
Dzyaloshinski-Morya (DM) interactions in the cluster 
Mni2. In our opinion, this can be attributed to the low 
symmetry of the cluster. Indeed, the strength of the DM- 
interaction is governed to a large extent by_asymmetry of 
crystal field acting on the interacting ionsEj. An instruc- 



the DM-interaction 
re surface of a magnet, 



tive example is provided in Ref. 2 
can emerge for ions located at t 
even though these interactions are prohibited for ions in 
the bulk of the magnet. In some sence the Mni2 molecule 
possesses " surface" everywhere, and the symmetry of the 
crystal field is rather low. 

The presence of the large Dzyaloshinsky-Morya term 
in the Hamiltonian provides a key to an explanation of 
the neutron scattering data. First, the DM-terms lead to 
the appearance of the two neutron peaks around 1.2 THz. 
If these terms are absent, there are two degenerate levels 
with S z = 9 around 1.2 THz. Among all the interactions 
we considered (see Sec. ||), only the DM-interaction can 
lift this degeneracy and provide a large splitting (about 
9 K), as observed in experiments. Similarly, according 
to our calculations, several peaks around 2 THz appear 
only due to DM-interactions. 

The origin of the peak at 0.3 THz has been completely 
explained in Ref. [l6|: it appears because of easy-axis 
anisotropy splitting the levels with different S z . Our re- 



G 



suits agree with this conclusion. 

An interesting feature in the neutron scattering data is 
the broad mode situated at 0.2 THz. No states of thk-ea.- 
ergy have been observed, e.g., in EPR-experimenta3'EJ'LJ. 
Our calculations also show no states with the energy 0.2 
THz (or, equivalently, about 10 K). In our opinion, this 
mode is caused by an interaction of Mni2 clusters with 
the dissipative environment. Due to this interaction, each 
level broadens (nonuniform broadening) , forming a qua- 
siband of finite width SE (see Refs. |ll|, ^ for details). 
The value of SE «2K can be estimated from the single- 
crystal hysteresis measurements^ Transitions between 
the two quasibands take place, so, along with the peak 
at 0.3 THz (14.4 K), a broad mode of inter-quasiband 
transitions appears. It happens when different states 
in the quasibands are populated, i.e. at the tempera- 
tures of order of SE w 2 K, which agrees with exper- 
iment. Energies of the inter-quasiband transitions are 
reduced by the value about 2 • SE ss 4 K, so the cor- 
responding neutron scattering mode is situated around 
E b = 14.4 K-2-6E « 10.4 K, or, equivalently, around 0.2 
THz, in agreement with experiment. At increasing tem- 
peratures, the occupancy of quasibands becomes more 
uniform, so the intensity of the broad mode increases 
along with the decrease in intensity of the peak at 14.4 
K. This behavior also agrees with experiment. However, 
this qualitative explanation can not be considered as suf- 
ficient, and a rigorous quantitative treatment is neces- 
sary. Such a treatment constitutes a separate physical 
problem to be investigated in the future. 

VI. DZYALOSHINSKY-MORYA INTERACTIONS 
AND THE RELAXATION OF MAGNETIZATION 

Strong Dzyaloshinsky-Morya interactions can play an 
important role in relaxation of magnetization in Mni2 
clusters. The possible importance of DM-interactions 
for the relaxational properties of Mni2 was already men- 
tioned in Ref. O, but no analysis was made. Below, 
the relation between the DM-interaction and relaxation 
is analyzed qualitatively. 

The relaxation time in M1112 is very long, about 2 
months at 2 K. At relatively high temperatures the mag- 
netization relaxation time r follows the Arrhenius law 
r = to exp (-AE/kT) with go = 2.4 • 10" 7 sec and the 
effective barrier AE = 61 KB. At low temperatures the 
situation changes: the experimental results presented in 
Refs. O, |l3| show that r saturates below the crossover 
temperature T c = 2 K, which is attributed to the onset of 
"ground state - to - ground state" tunneling. Such a high 
value of the crossover temperature poses certain prob- 
lems, since at present the nature of magnetic relaxation 
in Mni2 is unknown. Important infownation has been 
supplied by hysteresis measurementsol!9: it has been di- 
rectly shown that the relaxational mechanism in Mni2 
should allow for transitions between the states |<S Z ) and 
|5,±1>. 



The relaxation of magnetization has been extensively 
studiedt^EaEl within the single-spin model of Mn 12 . The 
Hamiltonian (||) of this model provides transitions with 
AS Z = ±4 only, and an external transversal field H x is 
necessary to allow the AS Z = ±1 transitions. Different 
sources of this field have beea considered: a hyperfine 
field induced by nuclear sninst3 and a dipole-dipole field 
induced by other clusterso. There are certain difficulties 
facing these interpretations. E.g., as reported in Refs. [ll], 
|3~i| , the relaxation rate increases on dilution of Mn 12 in 
solution, although the dipole-dipole interactions between 
the clusters decreases; this seems to be difficult to ex- 
plain with the relaxation mechanism based on interclus- 
ter dipole-dipole interactionscJ. The other mechanism, 
based on the hyperfine fieldsES can not easily explain the 
high value of the crossover temperature. 

Strong Dzyaloshinsky-Morya interactions constitute 
another source of magnetic relaxation in Mni2. It ob- 
viously allows the AS Z = ±1 transitions. On solution, 
the dipole-dipole fluctuating fields between the clusters 
decrease, thus decreasing the decohering influence of the 
environment, so the relaxation rate increases. This agrees 
with experimental results!! Relatively large magnitude 
of DM-interactions can, in principle, explain the high 
value (2 K) of the crossover temperature. This shows 
that the quantitative study of DM-based relaxation in 
Mni2 is important. 

Nevertheless, we emphasize that, to our knowledge, 
current information is not sufficient for judgement in fa- 
vor of some single relaxation mechanism. At present, 
all of them can be considered as equally probable, and 
the possibility of a combination of different mechanism 
exists. 

Independent of the problem of relaxation in Mni2, DM- 
interactions present an interesting and important mech- 
anism for magnetic relaxation. It can be significant, 
for example, in nanosized particles, where large DM- 
interactions can arise due to reduced symmetry at the 
surfaces. The relaxation based on this interaction is 
rather unusual: it is related not to the potential bar- 
rier created by easy-axis anisotropy, but to the barrier 
created by isotropic exchange. The point is that the 
Dzyaloshinsky-Morya term i?DM couples the states with 
different values of total spin S, i.e. 

(S,S Z \H DM \S±1,S Z ±1}^0. (13) 

This feature differs drastically from standard consider- 
ation of relaxation in small particles, where only the 
anisotropy barrier is usually taken into consideration. 

Isotropic exchange interactions are usually rather 
strong, but it does not necessarily means that the DM- 
based relaxation is negligible. E.g., in spin-frustrated 
systems the height of the exchange barrier can be con- 
siderably reduced. Moreover, it can be significant in cer- 
tain non- frustrated systems. For the relaxation based 
on DM-interactions, the relaxation time should be gov- 
erned by the ratio D/A, where D is the absolute value of 
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Dzyaloshinsky-Morya vector and A is the exchange inte- 
gral. On the other hand, for "conventional" relaxation 
mechanism, when the anisotropy barrier is overcome, the 
ratio U/K governs the relaxation rate, where K is the 
anisotropy constant and U is the strength of the interac- 
tion between the spins and their dissipative environment 
(e.g., spin-phonon coupling constant). Situations with 
D/A^>U/K are not impossible. 

Therefore, quantitative investigations of relaxation 
mechanisms based on DM-interactions is of great interest 
and importance. 

VII. SUMMARY 

In the present work, we have performed an extensive 
study of spin excitations in Mni2, explicitly account- 
ing for its many-spin internal structure. The dimer- 
ized 8-spin model of the Mn 12 clusters^ has been used. 
Along with isotropic exchange coupling, various kinds 
of anisotropic relativistic interactions have been studied: 
anisotropic exchange coupling between the cluster ions, 
single-site anisotropics of easy-axis and in-plane type, 
and various kinds of Dzyaloshinsky-Morya (DM) interac- 
tions. Surprisingly, most of these interactions play only 
a minor role. 

As a result, we propose a basic many-spin Hamiltonian 
which includes isotropic exchange couplings, single-site 
anisotropics of easy-axis type and DM-interactions be- 
tween the cluster spins. Three possible sets of parame- 
ters are determined from the relevant experimental data. 
The results of our calculations reproduce satisfactorily 
various experimental results, such as positions of neu- 
tron scattering peaks, high-frequency EPR data and the 
experimental dependence of the magnetic susceptibility 
on temperature. 

In particular, our results suggest rather strong 
Dzyaloshinsky-Morya interactions are present in the 
Mni2 cluster. We have discussed qualitatively the possi- 
ble relation of these interactions to the unusual magnetic 
relaxational properties of Mni2. We emphasized that 
DM-interactions present interesting relaxation mecha- 
nisms worth further investigations. 

ACKNOWLEDGMENTS 

Authors would like to thank A. K. Zvezdin, B. Bar- 
bara, D. Garanin for many helpful discussions. This work 
was partially carried out at the Ames Laboratory, which 
is operated for the U. S. Department of Energy by Iowa 
State University under Contract No. W-7405-82 and was 
supported by the Director for Energy Research, Office of 
Basic Energy Sciences of the U. S. Department of Energy. 
This work was partially supported by Russian Founda- 
tion for Basic Research, grant 98-02-16219. 



L. Thomas, F. Lionti, R. Ballou, D. Gatteschi, R. Sessoli 
and B. Barbara, Nature 383, 145 (1996). 
F. Hartmann-Boutron, P. Politi and J. Villain, Int. J. Mod. 
Phys. 10, 2577 (1996). 

12 B. Barbara, W. Wernsdorfer, L. C. Sampaio, J. G. Park, C. 
Paulsen, M. A. Novak, R. Ferre, D. Mailly, R. Sessoli, A. 
Caneschi, K. Hasselbach, A. Benoit, L. Thomas, J. Magn. 
Magn. Mater. 140-144, 1825 (1995). 

13 C. Paulsen and G. Park, in Quantum Tunneling of Magne- 
tization (Ref. p. 189. 

14 A. Caneschi, D. Gatteschi, L. Pardi and R. Sessoli, in 
Perspectives on Coordination Chemistry, ed. by A. F. 
Williams, C. Florani and A. E. Merbach (VCH, Basel, 
1992). 

15 A. K. Zvezdin and A. I. Popov, Sov. Phys. JETP 82, 1140 
(1996). 

16 M. Hennion, L. Pardi, I. Mirebeau, E. Suard, R. Sessoli 
and D. Gatteschi, Phys. Rev. B 56, 8819 (1997). 

17 A. M. Gomes, M. A. Novak, R. Sessoli, A. Caneschi and 
D. Gatteschi, Phys. Rev. B 57, 5021 (1998). 

18 B. Barbara, L. Thomas, F. Lionti, A. Sulpice and A. 
Caneschi, J. Magn. Magn. Mater. 177-181, 1324 (1998). 

19 A. A. Mukhin, A. K. Zvezdin, V. V. Platonov, O. M. Tat- 
senko, A. Caneschi, D. Gatteschi and B. Barbara, in Pro- 
ceedings of ICM'97 Conference, Cairns, Australia, July 27 
- August 1, 1997 (unpublished). 

20 Although this symmetry is not exact, the deviation from 
the symmetry is small. Such approximate symmetries are 
widely used in physics: e.g., the crystal lattice of ferromag- 
netic iron is often considered as cubic, though actually it 
is tetragonally distorted. 

21 A. L. Barra, D. Gatteschi and R. Sessoli, Phys. Rev. B 56, 
8192 (1997). 



8 



22 K. Yosida, Theory of Magnetism (Springer- Verlag, Berlin, 
New York, 1996). 

23 S. Hill, J. A. A. J. Perenboom, N. S. Dalai, T. Hathaway, 
T. Stalcup and J. S. Brooks, Phys. Rev. Lett. 80, 2453 
(1998). 

24 F. Luis, J. Bartolome, J. F. Fernandez, J. Tejada, J. M. 
Hernandez, X. X. Zhang, and R. Ziolo, Phys. Rev. B 55, 
11 448 (1997). 

25 See e.g., R. M. White, Quantum Theory of Magnetism 
(Springer- Verlag, Berlin, New York, 1983). 

26 A. Crepieux and C. Lacroix, J. Magn. Magn. Mater. 182, 
341 (1998). 

27 V. V. Dobrovitski and A. K. Zvezdin, Europhys. Lett. 38, 
377 (1997); L. Gunther, Europhys. Lett. 39, 1(1997). 

28 D. A. Garanin and E. M. Chudnovsky, Phys. Rev. B 56, 
11 102 (1997). 

29 A. Fort, A. Rettori, J. Villain, D. Gatteschi and R. Sessoli, 
Phys. Rev. Lett. 80, 612 (1998). 

30 A. L. Burin, N. V. Prokof 'ev and P. C. E. Stamp Phys. 
Rev. Lett. 76, 3040 (1996). 

31 R. Sessoli, Mol. Cryst. Liq. Cryst. 274, 145 (1995). 



Low-energy 


57 K 


58.2 K 


55.2 K 


56.7 K 


peaks (1.2 THz) 


66 K 


66.0 K 


66.7 K 


67.0 K 




maybe, 




67.4 K 


67.9 K 




more 


76.6 K 




75.7 K 


High-energy 


90 K 


124.4 K 


98.3 K 


88.8 K 


peaks (2 THz) 


96 K 


124.9 K 


105.1 K 


110.9 K 




105 K 


126.4 K 


110.4 K 


122.1 K 




maybe, 


127.1 K 


122.1 K 






more 









TABLE II. Parameters a and (3 of the easy-axis anisotropy: 
comparison between experimental data and calculated results. 
Calculations have been performed for the three possible sets 
of the cluster parameters (A, B and C, see text). 





Experiment 


Set A 


Set B 


Set C 


a(K) 


-0.56 


-0.63 


-0.68 


-0.67 


(3 (mK) 


-1.11 


-0.7 


-0.45 


-0.49 



FIG. 1. Schematic plot of the Mni2 cluster. Small black 
circles represent Mn 4+ ions, large white circles — Mn 3+ ions. 
Different types of lines connecting the ions (solid, dashed, 
dotted and dash-dotted) correspond to different types of ex- 
change interactions (Ji, J 2, J3 and J 4). 



FIG. 2. A schematic plot of the 8-spin system representing 
the Mni2 cluster. White large circles represent large spins 
(5 = 2), and dark small squares represent small dimer spins 
(s = 1/2). 



FIG. 3. Dependence of the normalized cross-section vs. 
level energy (in K), calculated for the three sets of the cluster 
parameters (A, B and C, see text). The levels producing no- 
ticeable neutron peaks can be easily discriminated from the 
others. 



FIG. 4. Temperature dependence of the effective magnetic 
moment of the cluster /j, e g (in Bohr's magnetons). Results of 
calculations with the three sets of parameters are shown: the 
set A (solid line), the set B (dashed line) and the set C (dotted 
line). Large solid squares represent experimental data. The 
results of calculations with the sets A and B are very close to 
each other, and the corresponding curves merge on the figure. 



TABLE I. The positions of neutron peaks: comparison be- 
tween experimental data and calculated results. Calculations 
have been made for the three possible sets of the cluster pa- 
rameters (A, B and C, see text). The levels with normalized 
cross-section more than 0.05 and energy less than 130 K are 
included in the table. 

Experiment Set A Set B Set C 
14.4 K 14.4 K 14.4 K 14.4 K 
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Fig. 1 

M.I. Katsnelson "Many-spin interactions ..." 




Fig. 2 

M.I. Katsnelson "Many-spin interactions ..." 



Parameter set A: J=0, J'=105 K, K 2 =5.69 K, 
D = 1 .2 K, D =25 K. 
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Parameter set B: J=23.8 K, J'=79.2 K, K 2 =5.72 K, 
D=-10 K, D =22 K. 
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Parameter set C: J=41 .4 K, J'=69 K, K 2 =5.75 K, 
D=-10 K, D =20 K. 
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Fig. 4 
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